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ABSTRACT 



Stellar density and bar strength should affect the temperatures of the cool (T ~ 20 - 30 K) dust component in the inner regions of galaxies, 
which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences 
between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 //m. We 
attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, 
which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey 
with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk 
regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk- 
temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient 
Ci , in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 ± 0.03. The cool dust temperatures in the central 
Q,, ~ 3 kpc of nearby galaxies are 23 (±3)% hotter for morphological types earlier than Sc, and only 9 (±3)% hotter for later types. The temperature 
i ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar 
density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, 
.£3 ' especially when paired with a strong bar. 

j3 Key words. Galaxies: ISM - Infrared: galaxies - Submillimeter: galaxies - dust, extinction 
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£> ■ 1. Introduction nificant star-formation activity, the resulting intense interstellar 

r~^ ' . , .. ,. . , radiation field (Stevens et al. 2005 ) and the earlier Hubble type 

iy-. The infrared emission from galaxies contains roughly halt of (-n Pnf i n Pt a i ?003i 

rv-, the entire energy budget in the Universe (e.g., Hauser & Dwek IT ' ... ' . c , cm c - A . CAA , 

CXJ ' ■ ---,-. T , ,. £ r b . . . b , r Here we use the infrared SEDs, from 70 to 500 i/m (a range 

_j uOOl). In addition to providing information on the amount of , . , , ,, , , . . , , .. c ■ tU 

1 ' . _ , , , ,, ,. , . r , • • which should be dominated by emission from grains in ther- 

• ' attenuation suffered by the stellar light, the infrared emission , .,.. • ... ., ,. .. c , , S * i 

iA • • , , • • , • • , , mal equilibrium with the radiation field, e.g., Popescu et al. 

£j , provides clues to important physical quantities, such as the 2AAA Engelbracht et al . 200 8), of a local sample of 13 galax- 

O metal, dust, and cold gas content of galaxies (e.g., Draine et . ,, ? — . , , . . ., TT , , , 
. -,™~, r, . . ™™ T jt . . .• •, ies spanning the range of spiral galaxies in the Hubble se- 
al. 2007 Bernard et al. 2008). Infrared spectral energy distnbu- . f ■ ., . . e ., , /rp OA -, A V \ 
^^ ' • ,ot^ s • „ i — ,• • , t •„• quence, to derive the temperature of the cool (T~ 20 - 30 K) 
• • , tions (SEDs), especially those extending into the submillimeter 1 , , ,,, , , . , ,, \. , , , 

>. , , , , dust component of the central region and the disks separately, 

regime, can be used to measure the dust mass and temperature , . \. . ,- a ■ ^. j « u .• -4.1.4. 

• ?-h , . b ' . . _ , _„„„ _„„, _ . , _„„. and investigate differences in the dust heating in the two re- 

k> in galaxies (e.g., Dunne et al. 2000, 2001; Seaquist et al. 2004; • ~, ■• ... , ., oeA ac:A ■ CAA 

^C '.,,,. v . & ,L„_ „ . , „„„„ ,,,.„ M , ~ nnn , . gions. To achieve this goal, we use the 250, 350, and 500 um 

r r\ ■ Vlahakis et al. 2005. Draine et al. 2007; Willmeretal. 2009, Liu %, , , /CDTDC . , c & 1 A T>t, t t - T : dc • 

H onirili Herschel/SPIRE (Spectral and Photometric Imaging REceiver; 

, . . 1 J; "■ , , , Griffin et al. 2010) images of the galaxies combined with the 

These temperature measurements have shown that the dust c ., /./r T Tic/i/r uu jt • m. , * e c -. r>- 1 

\ _ , . , . 5p;fzer/MIPS (Multiband Imaging Photometer tor Spitzer; Rieke 

is warmer in the centers of galaxies than in the outskirts (e.g., ' , OAA/l , -, A ji-- a ■ n * n •., papc 

,„ ^ 1 i , n »ni r. j ■ i \. , RnTiTi l i . i PtTTrTTn t^ et al. 2004 70 and 160 urn images. Eventually we will use PACS 

Alton etal. 1998, Radovichetal. 2001; Meloetal. 2002; Dupac ,„, t , t „ . JT 6 . c , \ r, v . , . , 

, i..-., ~ — ', , , ' — : — ' ' ' . r . (Photodetector Array Camera and Spectrometer; Poghtsch et al. 

et al. 2003). Cool dust at roughly the same temperature in spi- OA1A , • • f ./• -. j u * ,,u t - e+u- ,.■ tU 

, .., — ., ,,,,,, , , , nC r> , 2010) imaging for this study, but at the time of this writing, the 

ral disks is detected globally at longer wavelengths (850 um and TT — 7 * •, i_, tl i • u r„i , 

„. , ° J , i , -JL , _ b , i \„»„ i i^»», i data were not yet available. These galaxies have little nuclear 

1.3 mm; Siebenmorgen et al. 1999, Dunne et aLgOOO 2001, .. ., ... .,., .., , t ..? , wrpinni . • 

v „ , , . , rrTTTrti , i . , ., — — r activity that might heat the dust, with only NGC 1097 having a 

Vlahakis et al. 2005), but there is also some evidence of a . , , ., • • ., r- : . , . , , 

1 " . , , , strongly active nucleus, so this is the first study to cleanly sep- 

warmer temperature component associated with the central re- , ,/ , . , , ,. c . , , .. , . . 

r , r , ..... arate the tar-infrared properties of central and disk regions in a 

gions. Warmer dust temperatures tend to be associated with sig- , c . . • 

b ^_ b sample of normal galaxies. 

* Herschel is an ES A space observatory with science instruments pro- Until recently, little work has been done to dissect the dust 

vided by the European-led Principal Investigator consortia and with im- emission in galaxies into sub-galactic components, owing to the 

portant participation from NASA. general paucity of infrared images with the required angular res- 

** Hubble Fellow olution and to poor long-wavelength sensitivity. Some recent 
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work includes a study of the galaxy pair NGC1512/1510 (a tar- 
get also discussed in this paper), which finds that the dust tem- 
perature in the central region of NGC 1 5 1 2 is slightly higher than 
in the disk and that there is a significantly higher fraction of 
warm dust, in agreement with the center of NGC1512 being a 
starburst (Liu et al. l2010t . Work by Munoz-Mateos et al. ( 120091 1 
examines radial trends in dust properties in a number of nearby 
galaxies. 

The Herschel Space Telescope promises to yield a break- 
through in the study of subgalactic components in galaxies. This 
paper is the first investigation to leverage the longest wave- 
length Herschel data available for the KINGFISH (key insights 
on nearby galaxies: a far-infrared survey with Herschel; this 
program is largely derived from SINGS, the Spitzer infrared 
nearby galaxy survey by Kennicutt et al. l2003b sample of nearby 
galaxies, which will eventually total 61. Companion papers from 
the science demonstration phase for this program showcase the 
shorter wavelength imaging (Sandstrom et al. 2010) and the 
spectroscopic data (Beirao et al. 1201 Oi l. 

Here we present new SPIRE images acquired during the 
first few months of Herschel operations, in the context of the 
KINGFISH Open-Time Key Project. We divide each galaxy into 
two spatially-resolved zones: the circumnuclear region and the 
surrounding disk. Then we compare central temperatures with 
those for the disk. 



2. Observations and data reduction 

We observed 33 nearby galaxies with the SPIRE instrument on 
Herschel (Pilbratt et al. 2010) in the scan map mode as part of 
the KINGFISH program over the period of November 2009 to 
January 2010. They were observed in the scan mode out to 1.5 
optical radii, to depths of (3.2, 2.5, and 2.9) mJy/beam at (250, 
350, and 500) pm. They were reduced using the standard cali- 
bration products and algorithms available in HIPE (the Herschel 
Interactive Processing Environment; Ott 2010), version 1.2.5 or 
2.0.0 (whichever version was the latest available when the target 
was observed, since the SPIRE calibration is not sensitive to this 
range of HIPE versions). We modified the offset-subtraction al- 
gorithm in the pipeline, by requiring it to measure the offset only 
outside bright objects in the field (usually the galaxy that we 
targeted) — this procedure reduced a small negative offset (6% 
in the first galaxy we observed, NGC 4559) in the background 
in the same rows and columns as the source to an undetectable 
level. Otherwise, the data reduction was performed using the de- 
fault pipeline. 



3. Analysis 

Of the 33 galaxies observed, we selected 13 which had both a 
large extent (i.e., greater than several 160/500 //m beam diame- 
ters, or ~ 2') and a bright disk (brighter than ~ 0.5 Jy at 500 pm, 
excluding the nucleus). These galaxies are listed in TableQ] For 
the wavelengths with significantly smaller beams (i.e., 70, 250, 
and 350 //m) than the largest used here, we convolved the data 
to the SPIRE 500 pm resolution using the kernels described by 
Gordon et al. (2008) and updated by those authors for Herschel. 
For each galaxy, we computed masks that isolated the disk and 
central regions. The outer mask is an ellipse sized to encompass 
the disk where it contrasts strongly with the background, at 1 .2, 
6.6, 2.7, 1.4, and 0.8 MJy/sr at 70, 160, 250, 350, and 500 pm 
(where the SPIRE data have been converted to MJy/sr as de- 
scribed below), respectively. The inner mask is a circle centered 



Table 1. Galaxy sample and computed temperatures 



Name 


Hubble 


D" 


Spectral 6 


Tdis^ 


T c 
1 center 




type 


(Mpc) 


type 


(K) 


(K) 


NGC0628 


SAc 


7.3 


HII? 


22.7±0.4 


24.2±0.5 


NGC0925 


SABd 


9.0 


HII 


22.3±0.4 


23.9±0.5 


NGC 1097 


SBb 


19.1 


Syl 


23.4±0.5 


28.7±0.7 


NGC 1291 


SBO/a 


10.4 




18.3+0.3 


27.6±0.6 


NGC1512 


SBa 


14.4 


HII 


21.3±0.4 


27.3±0.6 


NGC3621 


SAd 


6.9 




21.9±0.4 


25.6±0.6 


NGC4559 


SABcd 


8.5 


HII 


22.7±0.4 


24.2±0.5 


NGC4594 


SAa 


9.4 


Syl.9 


21.3±0.4 


22.9±0.4 


NGC4631 


SBd 


7.6 


HII 


24.5±0.5 


27.7±0.6 


NGC5055 


SAbc 


10.2 


Transition2 


22.0±0.4 


25.3±0.5 


NGC5457 


SABcd 


7.1 


HII 


23.4±0.5 


23.9±0.5 


NGC6946 


SABcd 


6.8 


ffll 


24.5±0.5 


26.6±0.6 


NGC7331 


SAb 


14.9 


Transition2 


23.2±0.4 


26.2±0.6 



" As listed in NED: The NASA/IPAC Extragalactic Database (NED) 
is operated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and Space 
Administration. 

* from NED or Ho et al. [19971 

c includes uncertainty in the photometric calibration 



on the galaxy peak, with a minimum size chosen to be as small 
as possible (i.e., set by the MIPS 160 //m beam) and with a radius 
inversely proportional to the distance to the galaxy. Thus, a sim- 
ilar physical region of ~ 3 kpc diameter (in which one would 
expect to find structures like a bulge, bar, ring, and/or inner 
disk, e.g., see Erwin & Sparke 120021 Athanassoula & Martinet 
1980 and Knapen 2005) is sampled in the central region of each 
galaxy. Our results are not sensitive to a modest change in the 
aperture size, for example, similar results are achieved with a 
fixed central aperture size of 40" for each galaxy. A sample mask 
is shown in Fig.Q] 

We measured flux densities in these apertures for each of 
the 3 #ersc/z<?//SPIRE bands, as well as the Spitzer/MTPS 70 and 
160 pm bands. The Herschel/SPIRE data were calibrated in units 
of Jy/beam, which we converted to Jy by assuming beam sizes 
of 371, 720, and 1543 square arcseconds in the 250, 350, and 
500 pm bands, respectively (Swinyard et al. 1201 Oi l. We assign 
uncertainties as the quadratic sum of the data-dependent scatter 
and (by far the larger term) the uncertainty in the absolute cal- 
ibration, taken to be 5% and 12% for Spitzer/M TPS at 70 and 
160 pm (Gordon et al. 2007 and Stansberry et al. 120071 respec- 
tively) and 16% (the quadratic sum of the 15% absolute calibra- 
tion uncertainty from Swinyard et al. 2010 and a 2% uncertainty 
in the size of the beam) for SPIRE data. 

To each set of data, we fit a blackbody with a frequency- 
dependent emissivity of 1.5. Our conclusions do not change if 
values of 1 or 2 are used, although the temperatures increase 
or decrease, respectively. The uncertainties in the temperatures 
were computed via a Monte Carlo simulation, in which we per- 
formed 10000 trials where the photometric measurements were 
allowed to vary in a normal distribution with a standard deviation 
indicated by the photometric uncertainty. An example of fits and 
data points is shown in Fig. [2] where we can see that the 250 pm 
flux density is underpredicted by this simple model. This hap- 
pens frequently in our dataset, and may be a sign that multiple 
temperature components are present in the regions we measured. 
The values of the Herschel/SPIRE uncertainty indicated by the 
uncertainty in the absolute calibration were allowed to vary to- 
gether; i.e., the calibration uncertainty was assumed to be corre- 
lated among the three bands. In contrast, we allowed the MIPS 
uncertainties to vary independently; while not being strictly in- 
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Fig. 1. Example of the masking approach used, shown here at 250 fim on a cropped image of NGC 5055. The image on the left 
shows the full, reduced image, while to the image on the right we have applied a mask (the blank, white regions) to all data but the 
disk. Both images have been scaled logarithmically, from -0.01 to 5 Jy/beam. The contour levels are scaled logarithmically from 0.1 
to 5 Jy/beam. Each panel is 6' across. 
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Fig. 2. Far-infrared SEDs of the disk and central region of 
NGC 1512, each fit by a blackbody function modified by a /T 1 5 
emissivity. 



Fig. 3. Ratio of central dust temperature to disk dust temperature 
as a function of morphological type. The dotted line is drawn at 
a ratio of 1 (i.e., at equal temperatures) to guide the eye. 



dependent (see Stansberry et al. 2007), the 160 /jm calibration 
observations, in particular, are dominated by observational scat- 
ter and can be treated independently of the 70 //m band. The 
temperature of each component is taken to be the temperature of 
the best-fit model, with an uncertainty determined by the Monte 
Carlo simulation. 



4. Results 

For each galaxy, we computed a ratio of the central to disk tem- 
perature (from data at wavelengths longer than 70 fim) of the 
cool dust component. The average ratio of central-to-disk tem- 
peratures in this sample is 1.15+ 0.03. We have plotted this ratio 
against morphological type in Figure [3] where the average ra- 
tio is 1.23 + 0.03 for types earlier than Sc and 1.09 + 0.03 for 
later types. The trend is poorly described by a line, because it 
has a correlation coefficient of 0.66. We plot the ratio against bar 
strength (as defined by the morphological type) in Fig. [4] where 
the average ratio is 1.29 + 0.04 for strong bars and 1.09 + 0.03 
for weak bars. 

A significant (3<x) correlation between temperature ratio and 
Hubble type is observed in Fig. [3] although the correlation is 
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Fig. 4. Ratio of central dust temperature to disk 
ture as a function of bar strength, as indicated by 
type. The dotted line is drawn at a ratio of 1 (i.e 
peratures) to guide the eye. 
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only marginal, at 2cr, if the most significant point (NGC 1291) is 
removed from the correlation. This trend suggests that the dom- 
inant stellar bulges in early-type spirals are able to heat the cool 
dust in the ISM to greater temperatures than their late-type coun- 
terparts. The effect is most pronounced in the earliest type galaxy 
in our sample, NGC 1291, in which the dust in the central region 
is almost 50% warmer than the dust in the disk. In comparing 
this to the spectral types in Table [1] we have found no obvious 
correlation between the temperature ratio and the presence of a 
central nonthermal source in our galaxies. 

Bars may play an important role in driving the correlation of 
temperature ratio with morphology. As shown in Fig.|4j a strong 
bar is correlated with a high center/disk temperature ratio, al- 
though the scatter in the strong-bar ratios is large. This could 
occur through bar-induced starbursts, since mid-infrared colors 
are redder in barred galaxies, implying warmer dust tempera- 
tures (Huang et al. 1996, Roussel et al. 120011 . Barred galaxies 
also tend to have higher star-formation rates (SFRs) than un- 
barred galaxies, possibly because bars drive gas into the inner 
regions where it can pile up in the inner dynamical resonances 
associated with rings (Hawarden et al. 1986 ). The bar-driven gas 
transport has been shown in CO observations (e.g., Sakamoto 
et al. I1999I I. Stars will form quite readily in such conditions, 
where the gas transported by the bar fuels star formation (Erwin 
& Sparke 2002, Boker et al. 2008; Comeron et al. 120081 . In fact, 
the frequency of inner rings in early-type barred spirals (Hunt 
& Malkan [T999l Erwin & Sparke I2U021 could conspire to pro- 
duce the effect seen in our sample. Similarly, bars are more fre- 
quently found in bulge-dominated galaxies (e.g., Masters et al. 
1201 Oi l, which would also contribute to the trends we see. The 
concentrations of CO gas in the centers have been observed in 
early-type spiral galaxies with intense star formation (Koda et al. 
2005). In addition to a dependence on stellar density, the higher 
temperature ratios in the early types shown here could therefore 
be associated with bar-induced star formation. This is consistent 
with the result we get if the 70 yum fluxes are not included in the 
modified blackbody fits. Without them, the correlation between 
temperature ratio and either Hubble type or bar strength disap- 
pears; i.e., measurements at wavelengths shorter than the peak 
emission are required to determine the dust temperature accu- 
rately. 

5. Conclusion 

We used far-infrared data from Spitzer and submillimeter data 
from Herschel to compute separate SEDs for the center and disk 
regions of 13 nearby galaxies observed as part of the KTNGFISH 
program. We fit those SEDs (at wavelengths longer than 70 //m) 
with blackbody functions (modified by a frequency-dependent 
emissivity) to compute temperatures. On average, the cool dust 
temperature of the central component is 15 + 3% hotter than the 
disk. We find that the central temperature is higher than the disk 
by 20% to 50% in galaxies of type SO to Sb, but only 9% higher 
in later types. This ratio is also higher (at 1 .29 + 0.04) in strongly 
barred galaxies than in weakly barred galaxies (at 1.09 + 0.03). 
The data therefore indicate that the large (or "classical") 
grains that dominate the far-infrared and submillimeter emission 
are warmer in the centers of those galaxies with a substantial 
bulge and/or a strong bar. This may simply be caused by the 
higher density of the radiation field in the centers of early-type 
spirals, enhanced star formation due to the bar, or some combi- 
nation of the two. A cleaner separation of morphological compo- 
nents (perhaps with larger samples and/or less distant galaxies) 
and a more thorough assessment of the density of starlight and 



star formation activity, plus an evaluation of the impact of central 
nonthermal sources, may help separate these effects. 

The analysis presented here illustrates the power of Herschel 
observations in characterizating the spatially resolved distribu- 
tion of dust in nearby galaxies. This power will grow with the 
use of better-resolved far-infrared SEDs as measured by PACS 
(Poglitsch et al. l20101 l. which will let us measure smaller and/or 
more distant galaxies and determine radial trends of dust tem- 
perature. 
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